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Abstract 

1/-^ . We present the complete 2-loop renormalisation group equations of the superpotential 

^!^ . parameters for the supersymmetric standard model including the full set of -R-parity 

O . violating couplings. We use these equations to do a study of (a) gauge coupling unifica- 

Q>^ . tion, (b) bottom-tau unification, (c) the fixed-point structure of the top quark Yukawa 

coupling, and (d) two-loop bounds from perturbative unification. For large values of 
the R-parity violating coupling, the value of as{Mz) predicted from unification can be 



^~^- reduced by 5% with respect to the R-parity conserving case, bringing it to within 2a of 

r^ • the observed value. Bottom-tau Yukawa unification becomes potentially valid for any 

value of tan P r^ 2 — 50. The prediction of the top Yukawa coupling from the low tan /3, 

^ ■ infra-red quasi fixed point can be lowered by up to 10%, raising tan/3 up to a maxi- 

j^ ■ mum of 5 and relaxing experimental constraints upon the quasi-fixed point scenario. For 

heavy scalar fermion masses 0{1 TeV) the limits on the higher family AL ^ operators 

from perturbative unification are competitive with the indirect laboratory bounds. We 

calculate the dependence of these bounds upon tan/3. 



a) b.c.allanach@damtp.caiii.ac.uk 

b) adedes(9v2. rl.ac.uk 

c) dreiner@v2.rl.ac.uk 



1 Introduction 

The first grand unified theory (GUT), of the electroweak and the strong interactions was 
non-supersymmetric and the unification scale was of order Mx = 10^^ GeV ||l|. In order 
to connect the GUT predictions at Mx with observations at presently accessible energies, 
the renormalisation group evolution of the relevant parameters must be taken into account 
0. Postulating unification of the gauge couplings at a high scale leads after renormalisation 
to one low-energy prediction, e.g. the electroweak mixing angle sin^ 6w- In 1987, it was 
first found that in supersymmetry the prediction for sin^ 9\y is in agreement with the data, 
while in the Standard Model it is not [^, Q. This was spectacularly confirmed in 1990 
with the precise LEPl measurements of the gauge couplings constants |^, |^. This is the 
most compelling "experimental" indication for supersymmetry and has lead to a flourish of 
activity on unification and supersymmetry ^ §]• These studies focused on the minimal 
supersymmetric Standard Model (MSSM) which is minimal in particle content and couplings 
and conserves the discrete and multiplicative symmetry R-parityQ [^ 

We refer to this model as the i?p-MSSM, i.e. the Rp conserving MSSM. 

If we require a supersymmetric Standard Model which is only minimal in particle content 
the superpotential is modified to allow for additional R-parity violating (/?p) interactions 
which are given in full below in Eq.( |2.3| ). The superpotential includes terms which violate 
baryon number and separate terms violating lepton-number. In order to avoid rapid proton 
decay either baryon number or lepton number must be conserved but not necessarily both. We 
refer to a model which violates just one of these symmetries as an ^p-MSSM, [|^ . Symmetries 



which can achieve this are for example baryon parity and lepton parity |11, 12 

^^^(_)3B+25^ L,= (-)^+^^. (1.2) 

Thus, both the i?p-MSSM and the |?p-MSSM require a discrete symmetry beyond Gsm and 
are theoretically equally well motivated |T0|. 



1.1 R-parity Violation and Grand Unification 

Given the intense study of unification in the i?p-MSSM it is the purpose of this paper to study 
the gauge coupling unification in the /?p-MSSM. At first sight, it might seem unnatural to 
study unification within the ^p-MSSM, since |?p is not obtained in the simplest GUT models. 
In SU{5) for example, the dimension-four R-parity violating interactions are given by the 
operator 

^i^jXk, (1-3) 

^B: Baryon number, L: Lepton number, S: Spin. 



where i,j,k = 1,2,3 are generation indices, and ipi, Xk are the 5, and 10 representations of 
SU{5) respectively. The operator ( p..3|) contains all the cubic terms of Eq.( p.3|) , i.e. both the 
baryon- and lepton-number violating interactions. This leads to unacceptably rapid proton 
decay or unnaturally small couplings (~ 10^^^) and thus must not be present. In SO{10) 
and in Eq grand unification the dimension-four R-parity violating interactions are directly 
prohibited by gauge invariance. 

It seems R-parity violation and GUTs are incompatible. The reason is that any R-parity 
violating symmetry which is consistent with the bounds on proton decay, such as baryon 
parity and lepton parity in Eq. ([L^ ), assigns quarks and leptons different quantum numbers. 
But in GUTs quarks and leptons are in common multiplets and thus must have the same non- 
SU{5) quantum numbers. This contradiction is resolved once the GUT symmetry is broken, 
i.e. for energy scales below Mqut- Once the SU{5) symmetry is broken, R-parity violating 
terms can be generated which are consistent with proton decay. 

In general, we do not expect a GUT to be the final theory, it leaves many of the same 
questions unanswered as in the Standard Model. For example GUTs do not include gravity 
and therefore it should be an effective theory embedded in a more fundamental one, such as 
M-theory. This more fundamental theory will lead to a set of non-renormalisable operators 



at the GUT scale such as fl2 



^A^jXk^. (1.4) 

This operator is suppressed by a mass scale Mx ~ Mqut- Here, S is a scalar field in the 
adjoint representation of SU{b) and A; is a dimensionless coupling constant. ip^, ipj, Xk and S 
can be combined to SU{5) invariants in several ways. When S receives a non-zero vacuum 
expectation value SU{5) is broken and the operators in Eq.( |1.4D can generate a subset of the 
^p interactions in the superpotential (|2.3| ), which are consistent with bounds on proton decay 
p!3| . Models of this nature have been constructed for the gauge groups SU{5) []TT], |T^, |1^, |15[, 
SU{5) X [/(I) [|16l 0, and SO{10) ^ 



Below the SU{5) breaking-scale the operators (|1.4| ) are effectively dimension-four opera- 
tors. Their dimensionless coupling constant fc <S> /Mx will run, i.e. it will be renormalised 
and it will contribute to the running of the other couplings in the theory. Thus even though 
at first sight GUTs and R-parity violation seem inconsistent, this is not the case. Unless pro- 
hibited by a special symmetry, we expect to have R-parity violation via non-renormalisable 
operators in any GUT. At low-energy, this will manifest itself in (effective) tri-linear R-parity 
violating contributions to the superpotential. Above the GUT scale we will have an SU{5) 
symmetric theory with for example one unified gauge coupling constant. 

1.2 Unification and Fermion Masses 



One particular aspect of unification we will focus on below is the GUT prediction mb{Mi 



u) 



mT-{Mjj) which has been very successful [^^i Hi! HMl. We shall study the effect of R-parity 



violation on this prediction. If there are no non-renormahsable operators leading to effective 
fermion masses below Mqut, or if these operators are highly suppressed then we expect the 
Yukawa unification to still hold in the presence of R-parity violation. This can typically be 
achieved by a discrete symmetry but should be incorporated in a general theory of fermion 
masses (or Yukawa couplings) . If the non-renormalisable terms have the form 



w = hl{l:)x^x,K + Kf{T.)x^i^^h,, (1.5) 

the mass predictions are maintained. Here /?.„, h^ are the SU{5) 5 and 5 Higgs superfields, 
respectively, h^j, /i^j are general functions of the adjoint Higgs field. When SU{5) is broken 
and <S>7^ 0, the usual mass terms are generated. 

In the MSSM, if one requires the Yukawa couplings to unify this greatly reduces the 
allowed region of the (supersymmetric) parameters. In particular one obtains a strict relation 
between the running top mass mt{mt) and the ratio of the vacuum expectation values (vevs) 
of the two Higgs doublets, tan/9 []^, ^. Given the observed top quark mass |^ this results 



in a prediction for tan/5 ~ 1 — 3 or tan/3 ~ 55. Does this prediction still hold when allowing 
for R-parity violation? In Ref. ^^, by allowing only the bi-linear lepton number violating 
terms, it is shown that bottom-tau Yukawa unification can occur for any value of tan/3. The 
bi-linear term induces a tau-sneutrino vev, which introduces an additional parameter into the 
relation between A( and rrit, as compared to the MSSM. Bottom-tau Yukawa unification is 
then obtained by varying the stau vev, and therefore A^ (and hence tan/3). Here, we will 
focus on the effect of the tri-linear |?p terms upon the bottom-tau unification scenario. The 
third generation ^p-couplings enter the evolution of rrit, ruf,, and m,- at one loop and can thus 
have a large effect. Thus if we allow for ^p we expect the strict predictions of the MSSM 
to be modified. In Section ^ we shall analyse this effect and show that bottom-tau Yukawa 
unification becomes viable for any value of tan/?, each one corresponding to a particular value 
of an |?p coupling. 

There has been much work to predict the fermion masses at the weak scale from a simple 



symmetry structure at the unification scale ||T8| , [23|] . It is possible that the fermion mass 



structure is determined by a broken symmetry |2^ where only the top-quark Yukawa coupling 
is allowed by the symmetry at tree-level. Its value is put in by hand and is presumably of 
order one. The other couplings are then determined dynamically through the symmetry 
breaking model. Given such a model, we would then still require a prediction for the top- 
quark Yukawa coupling. An intriguing possibility is that this Yukawa coupling is given by 



an infra-red (quasi) fixed point ||2^ . The low-energy value then depends only very weakly on 
the high-energy initial value; the exact opposite of a fine-tuning problem. In supersymmetric 
GUTs with bottom-tau unification one typically requires large values of A^ ~ 1 close to the 
IR quasi fixed-point. This has been studied in detail in Refs.[0, |2^, |l^, ^ |2^. We investigate 



the effect of the |?p-couplings on the fixed point in Section |^. Similar to the case of bottom- 
tau unification in section E we find fixed-point structures for the top Yukawa coupling for 



any value of tan/3, although the focussing behaviour can be weakened depending upon the 
particular coupling introduced. 

1.3 Present Status 

There have been several previous studies of the renormalisation group equations (RGEs) of the 
^p- Yukawa couplings |^8|, |2^, ^ |3l|, ^, which have all been at the one-loop level. The main 



point of this paper is that we present the two-loop equations for the first tinie.0 In [pO| , |29[ 
the unitarity bounds on the couplings were determined at one- loop. These are still the best 
bounds on some of the baryon-number violating couplings. Below we update these bounds 
using the two-loop renormalisation group equations (RGEs). In |2^ the complete one-loop 
RGEs for the dimensionless couphngs were first presented and the fixed point structure was 
studied. We differ slightly in philosophy by also considering the Yukawa unification scenario 



as discussed above and considering the fixed point structure at two- loop in the RGEs. In |32] 
the full one-loop RGEs including the soft breaking terms were presented. These were used to 
study the bounds from flavour changing neutral currents. We do not here consider the RGEs 
for the soft breaking terms and this work is thus complimentary to ours. Several models have 



also been constructed implementing one- loop equations including the soft-terms |31|. Since 
our results are mainly model independent we do not comment on this work here. 

The most important effect which enters at two-loop is that the running of the gauge 
couplings now depends on the |?p-couplings. One might expect this effect to be small. But 
for higher generations the bounds on the /?p-couplings are weak and the couplings can be 
of order the electromagnetic coupling (e ~ 0.30) or more. In addition, most bounds are 
presented for scalar fermion masses of 100 GeV and become weaker for higher masses. At 
present the best la empirical bounds for the highest generation couplings and for relevant 
scalar fermion masses of 1 TeV are p!0| , |3i | 



A323 < 0.6 A'333 < 2.6 A"323 < 0.43*, (1.6) 

where the asterisk indicates the bound for a 100 GeV mass, and does not have a simple 



analytic description of the scaling with mass. At 1.5 TeV the bound on A323 ||3^ is almost 
identical to the perturbative limit obtained below in Section ^ The bound on A333 |^3 at 
1 TeV, is obtained by scaling and as such is meaningless since perturbation theory breaks 
down below Mu for such large values. The appropriate bound is thus the perturbative limit, 
which we obtain in Section |^. A mass-independent bound on A323 was found in ref. [^ from 



requiring perturbativity up to the scale Mu. However, we show below that the bound from 
perturbative unification is dependent upon tan/3, and we calculate this dependence. We shall 
thus explore all three couplings to the perturbative limit. 

^The two-loop equations have been presented before in |33]. This work contained a sign error in the 
RGEs as pointed out by the authors of [Q and remained unpubhshcd since one author left the field with the 
computer program. 



The outline of the paper is as follows. In Section 2 we present the full two-loop renormali- 
sation group equations for the gauge couplings and the superpotential parameters. In Section 
3 we present the specific RGEs assuming there is only one _^p-operator present at a time. In 
Section 4 we outline the procedure for our numerical analysis. In Section 5 we use the specific 
equations to study the effects of R-parity violation on the unification of the gauge couplings. 
We focus on the main predictions of unification: the unification scale, the value of the gauge 
coupling at the unification scale and the value of the strong coupling at Mz- In Section 6 
we study the effects of the third generation _^p-couplings on h-r unification. In Section 7 we 
study the Landau poles and the fixed points of the top- and the ^p-Yukawa couplings. In 
Section 8 we present our conclusions. 



2 Renormalisation Group Equations 

The chiral superfields of the i?p-MSSM and the ^p-MSSM have the following Gsm = SU{3)c x 
SU{2)l X U{1)y quantum numbers 



L: 



^-'-5'^) 2''' 



E: (1,1,1), Q: (3,2,1), f/ ; (3, i 2) 

D 3 



D: (3,1,-i), H,: (1,2,-^), H^ : (1,2,1). 



In the following we shall apply the work of Martin and Vaughn (MV) [^ to the general 
^p-MSSM superpotential. For the generic superpotential, W , we closely follow their notation 

W = F*''^'^<^*0,0.0,/6, (2.2) 

where (t)p,fj,& denote any chiral superfield of the model. The indices p, cr, 5 run over all gauge 
and flavour components. The ^p-MSSM superpotential is then given by 



W 



+eab 
+(^ab 



'y^E)ijL1H]^Ej + (Y£))ijQ'^^H-^^Djx + {Yu)ijQ1^H2Ujx 



a/^xb j 



{AEk)ijL1LjEk + {K]^k)ijL1Q^^ D^x 



+ -6.,.(A^o.fcf^r^p^ 



ra Tjb 



i T a Tjb 



fiH^H^ + K'L^H', 



(2.3) 



We denote an SU{3) colour index of the fundamental representation by x,y,z = 1, 2, 3. The 
SU{2)l fundamental representation indices are denoted by a,b,c = 1,2 and the generation 
indices by i,j, k = 1,2, 3. We have introduced the twelve 3x3 matrices 

Ye, Yd, Yu, A^^, A^., A^;,, (2.4) 

for all the Yukawa couplings. This implies the following conventions in the MV notation 



yi^QfS,, 


^ yLfSfc^Qf ^ yDkyL-Qb- ^ yQ^/L-Dky 






= yOf^^y^t = y5,,Q-Lr ^ (a^,),^.^.,^,^ 


(2.5) 


Y^tL'jEk 


= F^?^^^? = F^^^?^? = (A^0.,6a. = -{A^khe^,, 


(2.6) 


yUixDjyDk^ 


= YD.yU..D,^ = yD^yD^JJ.. ^ ^^^^ (A^O.fc = -exy. (A^.)., , 


(2.7) 



We denote the Gsm gauge couplings by g^, g2, gi- In Appendix A we have collected several 
useful group theoretical formulas pertaining to Gsm and the above field content. Here we 
mention that for f/(l)y we use the normalisation as in GUTs and thus use gi = Js/bgy- 
More details are given in the appendix. We define our notation for the Yukawa couplings via 
the superpotential including all ^p terms. 

We now in turn study the dimensionless couplings and then briefly also discuss the mass 
terms fi, Ki. We do not consider the soft-breaking terms here. 



2.1 Gauge Couplings 

The renormalisation group equations for the gauge couplings are 






IGtt 



2^2 



JlB^iSal- E fcTr(YtY.)+A^,ETr(AtA.j) 

6=1 x=uA,e \ i=l / 



:2.8) 



The coefficients !?„, Bab, and C^ have been given previously ||38l and for completeness we 
present them in the appendix. The ^p-effects on the running of the gauge couplings appear 
only at two-loop and are new. We obtain 



/ 12/5 14/5 9/5 \ 



A 



u,d,e 



V 




3 



6 

4 



1 
/ 



(2.9) 



This completes the equations for the running of the gauge coupling constants at two-loop. 



2.2 Yukawa Couplings 

In general the renormalisation group equations for the Yukawa couplings are given by |^ 



d 
~dt 



yijk _ yijp 



IGvr' 



^7^^^' + 



IGtt 



2^2 'P 



+ {k ^^ i) + {k ^-^ j] 



and the one- and two-loop anomalous dimensions are 



7. 



{i)i 



lY,p,Y^''^-2Sij:9lCa( 



(2.10) 



(2.11) 



7i 



(2)i 






+25|E^a 



glCa{i)Sa{R) + 2Y^glCa{i)CS) 

b 



?,glCS)C{Ga) 



(2.12) 



We have denoted by Ga{f) the quadratic Casimir of the representation / of the gauge group 
Ga- G{G) is an invariant of the adjoint representation of the gauge group G and Sa{R) is the 
second invariant of the representation R in the gauge group Ga- These quantities are defined 
in the appendix and their specific values are given there as well. We now first give the explicit 
version of Eq. (|2.10|) for the matrices ( p.4|) in terms of the anomalous dimensions, and then 
we present the explicit forms for ^yV \ and ^yV \ 



2.2.1 RG-Equations 

The RGEs for the Yukawa couphngs (including full family dependence) are given by 

[Y£;)ij = {YE)ikrEl + O^E)ijTHl - {^Ei)ki^Ll + C^E)kj'i^Ll^ (2-13) 



-^(Y^),, = {YD)^kro: + (YD)^JT^l-{Ao.)k^r'^l + {YD)k,r^l, (2.14) 



dt 
d 

di' 

|(Yf;),, = iYuUr^, + iYuy.,T^l + iYu)k,r'^„ (2.15) 

|(A^.)., = (A^O.^rf + (As.).;rg + {YEhkTi^, - (A^.),,ri; - (Y^),,r^^, (2.16) 



dt 

^(A^O.fc = (A^o.^r^r + (^^O/^r^; + (A^O.^r^:- (2.18 

At two-loop the anomalous dimensions are given by 



/. 167^2 '/' (IQtt 

2.2.2 Anomalous Dimensions 

The one-loop anomalous dimensions are given by 



r^; - jh-ff'' + 7T^7}!"'. (2.19) 



7i?^' 


= (YeY^e)^^ + iAE.Ai,),^ + 3(Az,,Al,,),, - Sfi^gl 


3 2^ 


) 


(2.20) 


■,r' 


= 2 (Yt^Y^)^.^ + Tr(A^. Aj,0 - 5/(^^?), 








(2.21) 


7S* 


= (Yz^yI,) ^.^ + (y^Y^)^.^ + (Al,,Az,0., - Sii^gl 


+ 


3 2 , 

2^2 + 


8 2, 

3^3)- 


(2.22) 


iT- 


= 2 (y^Yz,)^^. + 2Tr(Al,,A,5.) + 2(Af/,A[,,),. - Sji 


, 2 


2 8 

rgf + i^ 


■^3)), 


(2.23) 


4'!"' 


= 2 (Yt^Y,;)^^. + Tr(A,;. At,0 - SH^gf + ^^3^)), 








(2.24) 


1h, 


= Tr (3Y,,yL + Y^Y],) - {^gl + ^gl), 








(2.25) 


4r 


= 3Tr(YuYl)-i^^gl + lgl), 








(2.26) 


^r- 


= 7ffi^" = -3(A^,Yb),, - (A^,Y^),,. 








(2.27) 



Note that here, Hi2,L,Q represent the fields H^2j -^"5 Q^" where a is the index of the 
fundamental representation of SU{2) (i.e. no factors of tab are factored in). For the two-loop 
anomalous dimensions we write 

^P'' = h?'') . +(7^0 +(7^0 • (2-28) 

■" \ J^ / yukawa \ ■'* ' 9—y ^ J gauge 



These correspond respectively to the three terms of ( |2.12| ). These are given exphcitly below. 
The pure gauge two-loop anomalous dimensions are given by 



( 



iT 



The mixed gauge 

bT 



gauge 



0, 



gauge 

gauge \"± ^ gauge 

-Yukawa two-loop anomalous dimensions are given by 



(2)Uj 



bl 

(7£"" 



wsr 



b'S"' 



gauge 



gauge 



gauge 



gauge 



gauge 



rj/ ^4 -L^4 -'-'-''^4. o22 22 22\ 

^ii-g93 + -j92 + ^9i + 89392 + — c/3^1 + Y^929i) 

Tl / " 4 , -^UZ A oL 2 2\ 



bh 



84 856 4 128 2 2x 
'9^^ + 225^^ + 1^^^^^)' 



(-ri_-(7ro 



7g^- 



9-y 



9~-y 



9-y 



9~y 






91[[YdY^^ + 2YuY\j)^^ + (aL,A 



^Di /..J, 



J« 



,16 2 16 2n / » it 
(y^3 + Y^^i) (Ac7«A[;, 

+ (6^2 + \9l)[{YnY^^)^^ + Tr(A^. Aj^O], 



9-y 



9~y 



9~y 



(6^2' - \9l) (y^Y^)^^. + (^^1 - ^^?)Tr(A^. At,0, 
(16^3' - ^^?)Tr(Y^Yl,) + ^^,^Tr(Y^Yt,), 
{IQ9I + ^(^DTr (y^Y^) , 



(2)L>\* 



(7i,r') 



■|g; - 16g|) (K'„Yo)„ - ji (t.%,YE)„ 



'9-y ^ "" ^ 9-y 5 

The pure Yukawa two-loop anomalous dimensions are given by 

-(^rO,„._ = 2(Y,Yt,Y,Yy^.^+(Yt)^^(Y,),Tr(At,,A, 

+ 2 (A^^At,,)^.^ (y1.Y^),, + (A,,.ALO,.Tr (A^.A^, 

+ {YeY^)^^ Tr [YeY^ + SY^Yy 

-F (Y^).,(3Aj,,A^.YI, + A^,A^.Y^)^^ 

- (Y^),,(3A^.A1,,Y^ + A^.A^,Y^)J 

+ {K^j^kY eY^eKe^ + 3AjjfeAz)pAl)pA£;fc + AjjfcAi5pA^pA£;fe 



2.29) 

2.30) 

2.31) 

2.32) 

2.33) 
2.34) 
2.35) 

2.36) 
2.37) 
2.38) 

2.39) 
2.40) 
2.41) 
2.42) 
2.43) 



(2.44) 



7s, 



yukawa 



( (2 



(2)Qj 



yukawa 



iT- 



yukawa 






yukawa 



{i'«\ 



(2)Hi 



yukawa 



+ 3 ( A^feY/jY^A^fe + A^feY^Y^A^fe j . . + 3 f A^jk A^pA/jpA^^. 



J* 



= 2 (y1,Y^y1,Ys + y1,A^, Aj^^Ys + 3Y1,A^, A^, Y^ 

+ 2 {y^Ye) ^^ Tr (yI, Y^ + 3Yb Y^ 

— 2 (3Y^A£;i A^mY/) — Y^A^jj A^mY^; ) 

+ 2Tr [As. AJ,, (y^jYIj + A^. AJ,, + 3A^.A^, ' 



(2.45) 



2 (3Y^A£)mA^iY£; — Y^A^jmA^jY^; 



m] 



— 2 (^Y/pY^Y^Y^ 
+ 2{YuY\jYuYI 



n 



n 



j« 



Yz^Yl^ ,.Tr Y],Y,, + 3Yl5Yz5 



3(Yt/Y^) Tr(Y;,Y, 



J* 



(7^^t^ 



+ 2(a1,,A 



j^ll\£)m 



y 



Yl^Yz,)^ + (A^,At,, 



mi 



Im 



Tr(A|5,„A^! 



+ ( A^jmY^jY^A/Jm 1 + 3 ( A^jmAj).} A^jqA/Ji 

+ (Yb),,(3A^,Az,™Y1, + A5,,Ae™Y^ 



A.j-),n J^EI A^9 A/)n 



2 ( Y/jA^qAt/gYj) 



J* 



+ 2 (Y,,)^.„ (Yl, 



li 



TrfAl^^A^i 



(2.46) 



+ 



DJii 



^AI,A*o^Yd + A^^A^^Y^) .^ + {Yu)^i (Yi/),Jr (a^,A^. 



2 (yI^Yz^yI^Yz,)^ . + 2 (yI^Y^Y^Yz, 



+ 2 (yI^Yz,) . . Tr (y1,Yz, + 3Y1,Yz5) + 2 (y^AI,, Az^p Y^ 



J* 



+ 2Tr (A^, A^. (Yz^yI, + YuY\j) + A{,, Az,. A^Az,, 

+ Tr (6Az,.Al5,Az5. AJ,, + 2Az,,a1,,Az). AJ^, + 2Yz,y1;Az,. AJ 

+ 4 ( Az/m YJ, Yz) A[r™ ) . . + 4 ( A[;™ Az/p A[;p Ac/" 



j« 



(2.47) 



- 4 A 



"■tz-y.^l^c/" 



hi 



Tr (Az,.A^,) + (6Y1,A|;, Al,,Yz5 



ip 



+ 2 (y1,A^, A^pYz,)^^ + {qyIK^^^AdvY}, + 2Y|;a;,, Az^pY* 
+ 2(A^,Az 



3P 



Tr(A[,,A^fe 



2(Y^Yz; 



2 f yIYuY\jYu) . . + 2 (y^ Yz, Yl, Y^ 



4(AJ,,Az;. 



_l_ 9 /^v''" A^ A* Y 

+ 4Tr (Au.AI^Au.AI,) + 4Tr (aJ,,A^. yJ^Yz, 



6(Y^Yz; 



Tr {A^^rr^Aol 



= Tr (3Yz,y1,Yz,y1, + 9Y1,Yz,y1,Yz5 + 3Yz5Y1,Yz/YJ, 

+ Tr (3Yz;Y1;Az,, Al,, + Yz^Y^^Az,, A^,, + GYJ^Yz^aJ;, A 

+ 3Yz,y1,A^, Al,,) + (y1,Yz,) .. Tr (aJ,, Az,. 

+ 6(Yl,Yz,)^^Tr(At^,Az,; 



Tr (y^Y^) 



m 



(2.4^ 



(2.49) 
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\ ^ / yukawa \^ ~^ ^ ^ 

+ 3{yIjYu).Jt{A^u''Au.)' (2-50) 

+ (At,.Y^)^Jr(Al,,A^.) 

+ 3 I Y|;Y^ A^m Y^i — Y^Y^ A^m Y^; ) — 9 I A^jkYoY^Yn ) 

\ / im \ / ik 

+ 6 [a^^uYdAU^u^)^^ - 6 (Ad^Yz,),, Tr (a^. A^^ 



tA/_^ 'i ( s.* kT^ a* ^A_.^ _ _ /^,(2)^ 



- 3 A^,Yf/Y^YB ..-3 AI,,A^™A1,™.Y^ .. = - 7),^^^"M (2.51) 

\ / ij \ J ij \ ^ J yukawa 

This completes the renormahsation group equations for the Yukawa couphngs at two-loop. 
Before we discuss applications we briefly consider the renormahsation of the bilinear terms. 

2.3 Bi-Linear Terms 

Following the general equations given in MV the renormahsation group equations for the 
bilinear terms now including all R-parity violating effects are given by 

|a^ = /^{rS+rf^} + «:TjS (2.52) 

|«:^ = «:Tf^ + «:^r^; + /ir^^. (2.53) 

The anomalous dimensions at two-loop are given in the previous subsection. We see that 
even if initially Kj = a non-zero Ki will in general be generated through the RGEs via a 
non-zero fi. As noted in MV the bi-linear terms do not appear in the equations for the Yukawa 
couplings or the gauge couplings. They thus do not directly affect unification and we ignore 
them for the rest of this paper. 

2.4 Discussion 

The two-loop renormahsation group equations for the Yukawa couplings respect several sym- 
metries. If at some scale for example X'-j/^ = for all i, j, k then baryon parity, Bp, is conserved 
at this scale. There are no ^^p-couplings in the theory and thus in perturbation theory no 
^p-couplings are generated, i.e. the RGEs preserve X'-jf^ = at ah scales. Analogously, lepton 
parity, once imposed, is also preserved by the RGEs. If at some scale Xijk = X'-j^. = for all 
i,j, k and only one lepton flavour is violated, [e.g. A3 ;, 7^ in the mass basis of the charged 
leptons) then this is also true for all scales. If the neutrino masses are non-zero then this is 



no longer true ||53]. The electron mass matrix Y^; then contains off-diagonal entries which 



■^Ei -pLi 



generate off diagonal TJ. , T^^. via the RGEs. But the effects will be very small and can thus 
be neglected in most circumstances. However, if we assume only X'm 7^ at some scale then 
through the quark CKM-mixing the other terms X'^, will be generated. 
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Our results agree with MV for the MSSM Yukawa couphngs. We also agree with the 



one-loop Rp violating results |2^, |3§, ^ 



3 Specific RGEs 

3.1 Assumptions 

We now apply the two-loop RG-equations to the questions of unification. We shall assume 
as a first approximation that the ^p-couplings have a similar hierarchy to the SM Yukawa 
couplings and thus only consider one coupling at a time. The third generation couplings have 
the weakest bounds ( p.. 61 ) and can thus lead to the largest effects. We shall consider the three 
cases 

LLE : A323, LQD : X',,,, tJDD : A'3'23, (3.1) 

defined in the current basis of the charged fermions. We assume that in each case the re- 
spective operator decouples from the other /?p-operators whose couplings we set to zero. 
Rigourously this is not a consistent approach. However we can show that this is a very good 
approximation. 

The coupling (A£;3)23 violates L^ and can thus also generate non-zero (A£;i)i2 as well as 
{Aj^k)2j. We use (A£;3)23 as defined in a field basis in which Ye is diagonal, otherwise L^^r are 
not valid to protect the other (A^fc)jj from becoming non-zero. The full RGEs are directly 
determined from the equations in the previous section. The largest relevant terms at one-loop 
order are 

^(As0l2 ~ ^{KK){Ae^)23, (3.2) 

^(Ad3)23 ^ J^ihK)iAEs)23. (3.3) 

at IbTT^ 

Here he, h^-, hfy are the e, r, and b Yukawa couplings. We do not know the entries in Yd. Here 
we have assumed (Id) 33 ~ hb- Then in both cases the newly generated couplings are strongly 
suppressed and we can safely drop them. As an additional check, we have run the one loop 
RGEs in the above cases with (A£;3)23 non-zero at the unification scale. The initially zero 



M, 



V 



couplings are then generated at the 10~ level even for large values of {AE'i)22, 

The coupling {Kj^-i)^^ violates L^. and can in principle generate (A£,j)3j, (A^i)j3 7^ at 
the one-loop level. The relevant terms in the RGEs are 



-(A^.)3, ~ ^^-^(A^3)33 25,^{YlYDhk + h3[{YDYl)^, + {YuY^)^^] + 'iK{YD),k{^A) 



-(A^.)3, - Ye-TT 

1 Q 

-(A^Oai ~ Y^{Kh)iAD-^h3. (3.5) 






37(Ai52)32 ^ j^z^{hf,hb) {A 03)33, (3.6) 
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where j, k are not equal to 3 simultaneously. In order to estimate the contribution in the first 
RGE we would need full knowledge of the matrices Yd and Yjj, which does not exist to date. 
If we use the symmetric texture ansatze of [0 we indeed find the off-diagonal elements of Yd, 
Y^Yf) and Y^Yu highly suppressed. We thus feel justified in neglecting the newly generated 
couplings. In the second and third case we have made the same assumption as in ( p.3|) and 
we see that we can safely ignore the newly generated couplings. 

An initial non-zero (A{/3)23 can generate all the other baryon-number violating couplings. 
At one-loop and to leading order these are {Am)jf:, where (ijfc)e{(312), (313), (123), (223)}. 
The relevant terms in the RGEs are 



di 



(A 



WJjk 



i-K^ 



(A. 



t/3j32 



\YlYn)3A36j26ki + iY^YD)2A35jA3 + iY^Yuh^5,2h3] • (3.7) 



Under the previous assumptions this again leads to negligible new couplings. Thus in all 
three cases the decoupling is a good assumption. In line with this argument we assume the 
following form for the Higgs- Yukawa matrices 

YE = diag(0,0,A,), YD = diag(0,0,Afe), Y^ = diag(0,0. A*)- (3.8) 

So we make the approximation that the current basis is equal to the mass basis. 

3.2 Special Case Equations 

We now consider the specific case of the R-parity violation being through one dominant 
operator as in (p.l|). Thus we set the product of any two ^p couplings to be zero. We note 
the connection between the conventional notation and our matrix notation 



(A 



E'^ 32 



A 



323, 



(A 



D3J33 



A'. 



333' 



(A 



(73 J 23 



A 



323- 



(3.9) 



The gauge couplings have RGEs equivalent to those of the MSSM except for the following 
two-loop -^p contributions which we parameterise as 



A 



"• „ \ — ^M,-MSSM _ ^ Rp-MSSM 



dt 



dr" 



dr 



(3.10) 



where 






4 



(167r2)2 
9l 


(167r2)2 
9\ 



^^ 323 + ^-^ 333 



,'2 



6 A ooo + 2A 



2 
323 



(167r 



2A2 



L5 



- ^ 323 "T r; "^ 333 "T - "^ 



5 



5 



323 



(3.11) 
(3.12) 
(3.13) 
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The RGEs of the Yukawa couphngs are 

d 



dt 



-A 



323 



A323(rt + ri^ + r^^ +A.ri^^ 





= A'333 (rg + rg + rg) - A.r^^ 


d w, 


= A"3,3 (rg + rg + r^:) 


^A - 


- K {Til + rg + rg) + A323rS 


dt 


= A, (rg + ra + rg) - A'333r£ 


^A - 


= A, (rg + rt + rg) . 


The RGEs for the bi-hnear terms 


are: 


d 
Jt^ = 


/x (r^; + r^^) + /.^rg + /.grg 


d 

dt^' = 


-i(rii + rS) 


d 


/^r^^, + ^2 (rg + rg) 


d 


f^Tl\ + K, (rg + Til) ■ 



(3.14) 
(3.15) 
(3.16) 
(3.17) 
(3.18) 
(3.19) 

(3.20) 
(3.21) 
(3.22) 
(3.23) 



In the ROE for Ki the term proportional to /i is dropped because in our approximation T^l^ 
vanishes. The specific two-loop anomalous dimensions are: 



rg 



F^3 



+ 



1 



IGvr^ 
1 



2A^ + 2A 



323 



(167r2)2 
SAfcA 333 



6 2 

6 



(3.24) 



(A? + AL)(-5^? + %2) 
234 



6A1A 



2 \2 



4A323 



25 



■^1 



167r2 



At + A323 + 3A 



/2 
333 



10 



9l 



4A^ - 8A^A323 



-A 



(167r2)2 



15 



207 



-92 + 
4 ^^ 100 



^1 



9 

To 



gigi (3.25) 



+ A'L(16^3' - lal) + ^^?(A^ + AL) - 3A^ - 6A^Al23 

- 9A333A^ - 3AjA333 



'-'^323 



3A^Af, - 9A 333 



F^3 

^^1 



^ L2 






167r2 



3 A 333 A 



333/^6 



(IGtt 

1 
167r2 
1 



2^2 



A'333Afe(16(7| - -(7^) + 3A2A,A'333 - 9A'333AJ| - 9A 



'3 X 
333^b 



3A 333A. Ab 



'323 



10 



9l 



-A 



(3.26) 



(3.27) 



(167r2)2 



Y^2 + -^9X + Y5^2^1 + 5^lA323 



3A 



323 



3AtA323 



14 



pi/2 

'-Hi 









r^2 






P^2 






+ 



+ 



+ 



+ 



r 



Hi 
1 



167r2 
1 



167r 



2X1 



2-^323'^r + 



(167r 



2^2 



6 

— ( 



^323Xr-gi — 3A323Ar — 3A^A323 



(3.28) 



2A 



//2 



323 + ^^ 333 ~ 1 c^l ~ ods 



8 4 202 4 32 2 2 
9^3 + 225^^ + 45^^^^ 



(167r2)2 

16AfeA'333 - 8Af, - 2AfeAj - 8A 



(-^333 



(3.29) 



/4 
333 



2A 333AJ - 2A^A; 



U 



//4 
323 



4A 323At 



2A^A 333 



1 



167r2 
1 



^b + K + ^ 333 



(167r 

+ ^ 333) 
1 



2^2 



^ 33 

15 



30 



9l 



_ 2 _ _ 2 



(3.30) 



--^3 + —gt + 



— ^1 + 8^3'^2 + ^aWi + ^^2^? + ^^i(A6 + 2A? 



5A 



/4 
333 



lOA 333Afe 



\K 



//2 



5A{, — 5A4 — 2A 323A^ — A 



^/2 



333 ^V 



a! - 2A?A"' 



167r2 
1 



,//2 



(167r 



2 A 323 1 c^l 

_8 4 202 4^ 
9^^ 225^^ 



^3 



323 

(3.31) 



2A2 



^^«2 2 , X//2 /16 2 16 2n 4X//2 ,2 



^A 323 
1 



4A 323At 



167r2 

1 



3X1 + A' 



15 



10 



9! 



9I 



(3.32) 



(167r 



2^2 



-j92 + ^^^9t + l^9l9l + Xl{lQ9l-l9l) 



2\2 



:^lA^ 



3A^A323 - 3A^ 



3A^A 333 



9 A 333 Aft 



6A 323A5 



9A^ 



3AjA^ 



1 



167r2 
1 



3A? 



10 



9! 



9l 



(3.33) 



15 



207 



(167r2)2 

6 A 323 A J 



-92 + 

4 ^^ 100 



9i 



^9l9l + XK-/, + IQ9I) 



9Xt 



3A^Af, 



3 A 333 Aj 



167r2 
1 



2A 



- 8A 



(167r2)2 

//4 



323 + 2Aj - — fi'l - -^93 

8 4 856 4 128 2 '2,^ ,2/. 2 2 



(3.34) 



^?) + A" 



2 

323 



8 2^16 2 
-^9^ + y^3 



323 



4A 323Ab — 2AjA' 



333 



2A?Ag 



iX 



(3.35) 



4 Outline of the Numerical Analysis of RGEs 

In order to determine the scale of unification we numerically solve the renormalisation group 
equations. In the process, we re-derive unification predictions in the MSSM |0, ^7\, pO| , 



4T|-^^. The main aim of this paper is to isolate the new effects due to the SUSY flp 



part of the theory. When running the equations we must cross the mass thresholds of the 
supersymmetric particles. In order to determine these we must also run the two-loop RGEs of 
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Figure 1: Effect of non-zero |?p coupling A323 upon unification predictions. The value of tan/3 
input and the predictions for no ^p are shown in the figures, (a) Ri correspond to the ratio of 
the prediction of i in the Rp case to the predictions in the^p case, (b) Rb/r = \{Mu) / Xr{Mu) ■ 

the soft supersymmetry breaking terms. The l^p couphngs contribute to these RGEs, but this 
contribution has not yet been calculated at two- loop, although one- loop RGEs exist [^. In 
the future, when the full RGEs for the soft terms have been calculated, it will be interesting 
to include the /?p effects on the spectrum and thus on the thresholds. At this stage, we 
make the approximation of using the i?p-MSSM RGEs for the soft terms and the above ftp 
RGEs for the dimensionless couplings. We will also not consider (potentially large, but model 



dependent) GUT threshold corrections |^, ^ 



We add in turn one of the three /?p- Yukawa couplings (|3.1|). We run the full set of 
equations including the two-loop correction of the Yukawa couplings to the running of the 
gauge couplings in Eq. ( 



f). The boundary values of the running DR gauge couplings gi{Mz) 
and g2{Mz) can be determined in terms of the experimentally well known parameters, the 
Fermi constant Gp, the Z-boson mass Mz and the electromagnetic coupling a~^\j at Q^ = 0. 



a 



UQ' 



EM 



Gr p 

Mz 
= 0) 



1.16639 10"^ Gel/-^ 
91.187 GeV, 
137.036. 



(4.1) 
(4.2) 
(4.3) 



For a given set of pole masses rrit = 174 GeV, rrib = 4.9 GeV, rrir = 1.777 GeV we define the 
DR Yukawa couplings at Mz- Then we use 2-loop Renormalisation Group equations to run 
up to the scale Mu, where gi and g2 meetQ. 

■^The gauge boson self energies II^^, IIvi/iv are not modified by the R-parity violating couplings up to 



16 



1.1 



0.9 



0.8 



tan^=5 
a/M^)=0.128 



16 



M, =2.0x10 GeV 



a^=0.041 



(a) 




R 



R 



M, 



0.2 0.4 0.6 0.8 

^'333(111^) 



3.5 rr- 

3 - 

2.5 - 

2 - 

1.5 - 

1 - 



0.5 







(b) 



L^ 



..l-H-]- 



:r 



-/ 



0.2 0.4 0.6 0.8 

?l'333(m^) 



Figure 2: Effect of non-zero ^p coupling A333 upon unification predictions. The value of tan/3 
input and the predictions for no l^p are shown in the figures, (a) Ri correspond to the ratio of 
the prediction of i in the Rp case to the predictions in the^p case, (b) -Rb/,- = Xh{Mu) / \t{,Mu) ■ 

We have assumed universal boundary conditions at Mu for the soft breaking parameters 



Ao, Mo, Ml/2. The iteration procedure used to determine our results is described in [|J]. As a 



characteristic mean value of the soft breaking parameters at Mu we choose Aq = M1/2 = Mq = 
300 GeV. We solve the RGEs for different values of the ^p-coupling at rrit, starting from zero. 
The maximal value we consider is where the running coupling reaches the perturbative limit 
at the unification scale. For most of our discussion we take this to be A(M^) < -\/47r ^ 3.5. 
However, we also show some results where A(M[/) < 5.0. 

There has been a large interest in the literature ^ |20|, |1^, |18|, [19|, ^ in the restrictions on 
the unification scenario from bottom-tau unification. Requiring bottom-tau unification leads 
to a strict relation between the running top quark mass and tan/3. For the experimental value 
of rrit [0, tan/5 is predicted to be very close to 1.5 or around 55 [0. We are interested in 
how the effects of fip can relax this strict relation and allow a larger range of tan/5. As a 
model scenario, we consider tan/5 = 5 which is well away from the solutions in the MSSM. We 
show that this feature is general and that solutions may be obtained for any value of tan j3, 
provided one tunes the value of the ^p coupling. The ^p couplings change the fixed-point 
(and quasi-fixed point) structure of the top- Yukawa coupling's evolution. This is investigated 
in detail for low values of tan/5. 



1-loop level. We assume that the effect of those couplings in vertex and box corrections to the running weak 
mixing angle is negligible Wa. H. 
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Figure 3: Effect of non-zero ^p coupling A323 upon unification predictions. The value of tan/? 
input and the predictions for the no fi^ case are shown in the figures, (a) Ri correspond 
to the ratio of the prediction of i in the R^ case to the predictions in the /?p case, (b) 

Ri,lr = Xb{Mu)/Xr{Mu). 



5 Gauge Coupling Unification 

For A^p = and with the inputs and procedure defined above, we obtain 
a,{Mz) = 0.128, Mu = 2.0 lO^*^ GeV, au = 0.041. 



(5.1) 



In order to discuss the effects of the non-zero ^p Yukawa couplings we consider the unification 
parameters as functions of A^^^ evaluated at m^. 



and define the ratios 






(5.2) 



Ruui^np) 

RauiXftp) 



«3(M^,0) 



(5.3) 



auiO) 



When A^ = and tan/? = 5, the ratios in (|5.3|) are then all equal to 1 as can be seen on the 
left of Figs. |l]a, ^, and ^, respectively. Next we turn on the |?p-couplings. In Fig.s |l]b, ^, 
and ^ we can read off the value of the ^p-coupling at the unification scale as a function of 
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the coupling at rrit. The plots stop once the perturbative limits are reached. For the present 
numerical discussion we focus on tan/? = 5. X323{Mu) reaches its perturbative limit for a low 
scale value of A323(Tnt) = 0.93. It is worth pointing out that this is the same as the laboratory 
bound for slepton masses at 1.5 TeV! Thus although the laboratory bounds on the LLE 
operators are generally considered to be very strict; for heavy supersymmetric masses they 
are no stricter than the perturbative limit. At this point A323 has run off Fig. |l|b but it should 
be clear how it extrapolates. The perturbative limits for the other couplings are given by 

A'333(rrit) = 1.06 (5.4) 

A"323K) = 1.07 (5.5) 

The first limit ( p.4| ) is equivalent to the empirical 2 a limit for 700 GeV squark masses. The 
la empirical bound on A"323 is 0.43 for 100 GeV squark mass, and so the limit (|5.5| ) will be 
restrictive for somewhat higher masses. These limits are tan/3 dependent, and we leave a full 



discussion to section 1773 



In Fig.s |l|a, 0a, ^ja we show how the ratios ( |5.3| ) change as we turn on the ^p-couplings. 
For A323 7^ 0, as{Mz) and au are practically unchanged except very close to the perturbative 
limit. However, Mu is shifted upwards by up to 10(15)%, where the parenthesised value 
corresponds to choosing the perturbative limit to be A < 5.0. For A'333 the downward shift 
in as{Mz) is typically 1-2%. At the extreme perturbative limit A < 5.0 the maximum shift 
is is a decrease of 5% in as{Mz) giving a value of as{Mz) = 0.122. This corresponds to an 
agreement with the data as{Mz) = 0.119 ± 0.002 [^ to 1.5a, without using GUT threshold 
corrections. While the importance of this result is obvious, caution in its interpretation 
is required because the correction to as{Mz) is sensitive to where one places the limit of 
perturbative behevability. For example, if one chooses X{Mu) < 3.5, one only obtains a 3% 
decrease, still with significantly better agreement with the data than the Rp conserving case. 
au is decreased slightly at this point. However, Mu is decreased by up to 20%. This effect 
is significantly beyond the effect due to the top quark Yukawa coupling. For A"323 7^ 0, au 
remains practically unchanged, as now has an overall increase of up to about 3% at the 
perturbative limit corresponding to a value of as{Mz) = 0.131 in disagreement with the 
experimental value. Mu is raised by up to 20%. 

Thus we find au essentially unchanged by ^p-effects. Mu can change either way by up 



to 20%. If we compare this with other effects considered in Ref. |^ we find it of the same 
order as the uncertainty due to the top quark Yukawa coupling or the effects of possible non- 
renormalisable operators at beyond the GUT scale. The effect is much smaller than that due 
to GUT-scale threshold corrections or weak-scale supersymmetric threshold corrections. It is 
thus much too small an effect to accommodate string unification. The strong coupling can 
also change either way by up to 5%. A decrease is favoured by the data and is welcome in 
supersymmetric unification. The effect of the ^p-couplings on as{Mz) is of the same order as 
the effects due to the top-quark Yukawa coupling, GUT-scale threshold effects and high-scale 
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Figure 4: Bottom-tau Yukawa unification including the R-parity violating couplings (a) A333, 
(b) A323 for various tan/? = 1.6,5,26. The plots show i?^/^ = \{}vlij) j \t{Mu) as a function 
of the |?p couplings evaluated at m^. In (b), the tan/3 = 1.6 curve stops as \t{M\j) reaches 
its perturbative limit. 



non-renormalisable operators M . 



6 b-T Unification 



In order to study the unification of the bottom and r Yukawa couplings A;,, A^- at Mjj we 
define the ratio 

A5(Mt;,A^J 



RhiAMx 



u 



(6.1) 



A,(Mc;,A^J' 

For Xfi^ = 0, tan/3 = 5 we have 

Rb/r{Mu) = 0.78. (6.2) 

Thus including the top-quark effects but before turning on the ^p-coupling we are well away 
from the bottom-tau unification solution Ri,/t-{Mu) = 1. Recall that the uncertainties due to 
the bottom quark mass are small for small tan/?. Now we consider the corrections due to the 
|?p-couplings. The one-loop RGE for Rh/r(t) is given by 

Wn'^^^ = R,/r{t) [\l + 3A^ - 3A^ - 4A^,3 + ^\%, + 2A;13 + ^^? " f ^3^] • (6-3) 

The leading dependence of R^/r on A323 has a negative sign and as we see in the two-loop 
result shown in Fig. |I|b Rb/r drops significantly. Near the perturbative limit it drops by a 
factor of 2.5 and ^p becomes a dominant effect on the evolution of -Rb/r- This is important 
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for the range of tan/3 which leads to bottom-tau unification. In the MSSM R^It- is too large 



for tan/? ~ 1.5 or ~ 55 [^. Including a non-zero operator A323 strongly reduces /?&/,- and 



thus can lead to bottom-tau unification in this previous regime. 

For A"323 7^ or A'333 7^ there is an additional positive contribution in the evolution of 
Rb/T-{t). The full two-loop result shows a clear rise in RhfT-{t) as a function of A"323 in Fig. ^. 
The maximum increase at the perturbative limit is by 60%. For A"323(r?7,i) = 0.9 bottom- 
tau unification is restored! This is quite remarkable. Even though ^p-couplings are usually 
expected to lead to only small effects they can have a significant impact on our understanding 
of Yukawa-unification. Recall, that grand unification is possible in ^p-theories as discussed 
extensively in the introduction. From Eq. ( |6.3| ) it should be clear that for example for A'333 
we get an increase in Rh/T-{Mu) as well leading to further bottom-tau unification points. From 
Fig. II, we see that for A'333 (m^) = 0.75, we achieve Rb/r = 1 for tan/3 = 5. To investigate 
further, we plot Rb/r as a function of A'333 (rrii) for several other values of tan/? in Fig. ^. The 
figure illustrates that any value of tan/5 = 2 — 26 achieves bottom-tau Yukawa unification, 
provided A'333 (m^) is chosen correctly. Fig. Qd shows the equivalent plot for A"323 7^ with 
the same conclusions. 

7 Landau Poles and Fixed Points of Yukawa Couplings 



True fixed points and quasi fixed points were first considered in |^. Since then, many ap- 
plications have been discussed in the literature (see ref . fSO] and references therein) , including 
those in the Rp conserving MSSM. In the MSSM and at low tan/3, it is well known that if 
one neglects two-loop corrections as well as those from couplings smaller than g^, the RGE 
for Xt has an infra-red stable fixed point |^ corresponding to Xt/gs = a/7/18. This is too 
low to accommodate rrit = 175 ± 5 GeV. To be phenomenologically viable, A^ must be out of 
the domain of attraction of the true fixed-point. In practice, this implies that Xt{mt) must be 
nearer to its quasi-fixed point (QFP) limit of 1.1, where Xt{Mjj) is large (formally it diverges). 
Using mt{mt) = 165 ± 5 GeV as an input means that one can derive tan/? in this scenario 
through the relation 

. V2mt{mt) 

sm(3 = — — — — , (7.1) 

vXtimt) 



implying tan/? = 1.7 ± 0.2 at the QFP |^. This scenario is very attractive [p6| , Q because 



the values of many parameters in the infra-red regime are insensitive to their input values 
at Mu, giving higher predictivity and a tightly constrained phenomenology. The quasi-fixed 
Rp conserving MSSM is presently severely constrained [0, Q. A large |?p coupling changes 
the running significantly, and so in this section, we examine the running of the couplings as 
a function of renormalisation scale. It is convenient to split this analysis into three cases to 
focus upon: (1) infra-red stable fixed points of approximated RGEs, (2) a small amount of 
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^p near the QFP of At and (3) the case of two large Yukawa couphngs. We then examine the 
constraints on the ^p couphngs defined at nit coming from the requirement of perturbativity 
up to Mf/. 

7.1 Fixed Points 

First, we analyse the fixed points in the equations including the couplings A333, A323 one at a 
time. Neglecting A^, A,- (i.e. looking at the low tan/? limit) and two-loop terms, we reformulate 
Eqs. ( PTT^ , (p:TlD-(p:TBD as 

— ; — ^ = X' -2Xt--X" + — + 7.2 

dhigl 9 3 3 03 45^3 

d\nX 3 tti ao 4 ^^ ,_ , 

dmg^ 5 03 0:3 3 

dhigl 9 3 03 45 03 ^ ' 

d m 5f| 3 3 5 0:3 

where Xt = Xf/g^, X = X\2'i/9h ^' — -^'333 /S'l ^"^^ ^" — ^"12?,/ 9l- The stability of Yukawa 
couplings in supersymmetric theories has been considered in refs. [^. Considering only one 
non-zero l^p coupling at a time, we have two infra-red stable fixed points in the limit that we 
ignore the electroweak gauge couplings. The first is 

X = X' = 0: X, = \, X" = g (7.6) 

and the second is 

X = X" = : Xt = X' = -. (7.7) 



We have ignored A323 in this discussion because it does not exhibit fixed point behaviour 
itself due to the lack of renormalisation from the QCD interactions. The values of Xt in 
Eqs. (|7.6|),([7l^) are even lower than the _Rp-conserving MSSM value Xt = 7/18. They are 
experimentally excluded |5^ by the lower bound on \t{mt) coming from the requirement of 
'mtinrit) = 165 GeV. 

7.2 Large Xt{Mu), Small ^p Coupling 

We now discuss the case where Xt{Mu) > 1 and the ^p-couplings are in turn very small. The 
behaviour of the Rp conserving superpotential parameters will be similar to their behaviour 
in the Rp conserving MSSM. In this case, we can solve the one-loop RGEs for the ^p Yukawa 
couplings analytically: 

h23{fi) = \,2:,{Mu)5f5f'^\ (7.8) 

A'333(/.) = X'...{Mu)5T'5-^'h-^'"''5l'\ (7.9) 

A"323(/x) = \",,,{Mu)5f"'5,5-,'"''hl'\ (7.10) 
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Figure 5: (a) Effect of large R-parity violating couplings upon the quasi-fixed point of Aj. 
The ^p couplings increase toward the bottom of the plot. The lower limit of Xt{Tnt) is defined 
as the minimum value required to obtain mt{int) = 165 GeV. (b) Perturbative limits of 
large R-parity violating couplings as a function of tan/5. The upper curves correspond to the 
two-loop order calculation and the lower ones to a one-loop calculation. 

where ^1^2,3 = 9i,2,?,{lj) / 9i.^u) and 5t = Xt^jj) I \t{Mu) . We have neglected contributions from 
At-, Ab in Eqs. ( [7.8|) -( fnoD and so they are valid only at low tan/5 < 15. The one-loop analytic 
solutions for 5i,2,3 are equivalent to the Rp conserving MSSM ones 0]. The solution for 5t 
is also equivalent to its Rp counterpart, which has been solved analytically at one loop order 
including (73,(72 [HI- Eqs. ( l7r^ )- (|7.10| ) predict a constant ratio of Xfi^{Mu)/\^^{'mt). This 
corresponds to a straight line for the X^^{Mij) curves in Figs. [I|b,^,|^. For Xf^^imt) ~ 0.3 
this is a good approximation, after that the two-loop effects become important. We also obtain 
analytic solutions at low tan/3 for the bi- linear ^p terms when the |?p Yukawa couplings are 
small: 

n.{^^) = ^.{Mu)5-,"'''5-,''%"'5]'\ (7.11) 



7.3 Two large Yukawa couplings 



The one-loop RGEs including Af, A^, (72,3 have been solved in the MSSM |^. These RGEs 



have the same form in the case of a non-zero A333, when ignoring all other Yukawa couplings 



except At. The analytic solutions to this system are therefore contained in ref. [^ and are 
in terms of hypergeometric functions. 

We present here numerical solutions to the two-loo^) equations, which can also be applied 
to the cases of large A323, or A"323 with large Xtirrit). Using au = 0.041, Mu = 2 x 10^^ GeV 
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Figure 6: Running of Ai(/i) and (a) A323(/i), (b) X'^^^{fi). The circles are separated by two 
orders of magnitude of /i/ GeV. The flow is towards the region (a) A^ ~ 1, A"323 ~ 1 and (b) 
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1 for decreasing /i. 



and again switching off A^ and A,- (which is only valid for low tan/3 ~ 15), we calculate how 
\t{mt) is related to its input value at Mjy, i.e. we study the quasi-fixed point structure. SUSY 
threshold corrections are not included. 

In Fig. ^, the top solid line shows the quasi-fixed point structure of A( in the MSSM. 
For \t{Mjj) > 2.0 it is almost flat, i.e. \t{mt) = 1.1 becomes insensitive to Xt{Mu). The 
horizontal dashed line corresponds to the minimum \t{mt) required to produce a top quark 
mass which agrees with the data and for which sin/? < 1. It changes by ±0.03 within the 
empirical errors on mt. The A323 7^ curves are not plotted because they coincide with that of 
the MSSM. This can be understood from Eq. ( |3.19| ) as the coupling does not directly appear 
in the running of A(. This holds at two- loop due to our assumption of only one dominant l^p 
coupling at a time. 

When A'333 is switched on, the quasi-fixed point behaviour of A^ persists but its QFP 
value can be decreased as far as Xt{rnt) = 1.03 which corresponds to tan/5 = 2.1±q;p Any 
decrease in \t{mt) can increase tan/3 as extracted from Eq. ( [7.1D , modifying the QFP pre- 
dictions of superpartner masses, for example. In particular, higher tan (3 can allow for higher 
masses of the lightest CP-even Higgs, relaxing a severe constraint in the Rp conserving QFP 
scenario [^. For large A"323, the quasi- fixed behaviour of A^ is somewhat weakened as the 



non-zero gradient of the curves suggests, but Xtijrit) = 1 is possible as the QFP prediction^, 
corresponding to tan/5 = 3.0±q;7. 

Fig. 1^ only shows information on the quasi- fixed point structure of A^. One would like to 



*This lowering of Xt{mt) was first pointed out by Brahmachari and Roy in ref. pi 
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know if the ^p couplings also exhibit the QFP behaviour when they are large. Fig. ^ shows 
the running of both A((/i) and \'i2-ii.l^) with the renormalisation scale /x. A;,, At have been 
switched off in this calculation. The two-loop -^p-MSSM RGBs were run through 14 orders of 
magnitude (roughly corresponding to running from Mjj to nit), using au = 0.041. Although 
the figure shows both couplings running toward 1 in the infra-red regime, it is clear that it is 
difficult to make this statement quantitatively accurate. The same conclusion holds for large 
A333 7^ and At, as shown in Fig. ^d. A323 exhibits even less focussing behaviour because it is 
not directly affected by QCD interactions. 

7.3.1 Perturbative Limits 

In Fig. Qo, we show the limits from perturbativity upon the ^p-couplings defined at m^. We use 
a degenerate effective SUSY spectrum at rrit in this calculation but no finite SUSY threshold 
effects are included. We include A,-, A;, however, because they make a large difference to any 
Landau poles of Yukawa couplings at high tan/5. When we switch an fip coupling on, the 
curves in the figure map out what value of the coupling is required at rrit to produce a value 
of 5 for at least one of the Yukawa couplings at Mu. In practice, to good accuracy the point 
where a Yukawa coupling is 5 is very close to the Landau pole of that coupling. The deviation 
between the two curves shows the significant weakening effect of including two-loop terms in 
the RGEs: 5%, 12% and 10% for A323, A323 and A323 at high tan/3 respectively. The upper 
bound calculated in this way exhibits a strong tan/5 dependence due to A^, A,- contributing 
to the running at high tan f3. Our one-loop bounds in Fig ^ agree with the previous limits 
upon A323 from perturbativity bounds provided in by Brahmachari and Roy [^ ] to within 
3%. One-loop bounds upon other A"^^ were obtained by Goity and Sher ||36|1 . 

8 Conclusions 

We have argued that /?p is theoretically on equal footing with conserved Rp. Since it can 
be realized in grand unified theories it is relevant for unification. We then first determined 
the complete two-loop renormalisation group equations for the dimensionless couplings of the 
unbroken supersymmetric Standard Model. It is only at two-loop that Yukawa couplings 
affect the running of the gauge coupling constants. We then considered three models of 
/?p. We have added to the MSSM in turn the three Yukawa operators L3L2-E3, ^3(531)3, 
and U^D2D'i. We considered their effects on various aspects of the perturbative unification 
scenario. We have focused on qualitative effects. A detailed search for a preferred model is 
beyond the scope of this paper. We found several important effects. The unification scale is 
shifted by up to ±20%. This is comparable to some threshold effects but insufficient for string 



^Note that Brahmachari and Roy differ by a factor of two in their convention for the A" superpotential 
terms. 
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unification. as(Mz) can be changed at most by ±5%. The reduction which is favoured by 
the data is obtained close to the perturbative hmits of A'333 and A323. We have obtained the 
two-loop limit from perturbative unification for all three operators. For A323 it is equivalent 
to the laboratory bound for a slepton mass of 1.5 TeV and for A333,A323 is competitive for 
masses below 1 TeV. Two-loop limits from perturbativity are 5 — 12% weaker than the one- 
loop limit previously obtained. This is all quite remarkable. The /?p-couplings can have 
significant effects on the entire Yukawa unification picture. For bottom-tau unification we 
have found significant affects. For A323 7^ bottom-tau unification could be obtained for 
values of tan/5 < 1.5 were it not for the fact that the perturbative limit is reached. For 
A323, A333 7^ 0, we found new points of bottom-tau unification at tan/? = 5 — 26. Thus for 
tan/5 < 30, bottom-tau unification doesn't necessarily correspond to top IR quasi fixed-point 
structure, as is the case in the MSSM 0, |2^, |2^, |5^. Indeed, the quasi-fixed structure is 
changed resulting in lower values of the Aj(mj) prediction for A333, A323. This allows tan/? to 
increase, resulting in higher masses for the lightest CP-even Higgs, relaxing severe constraints 
on the QFP. 

Acknowledgements 

We thank Stefan Pokorski, Vernon Barger, Probir Roy and Marc Sher for helpful conversa- 
tions. This work was partially supported by PPARC. A.D acknowledges the financial support 
from the Marie Curie Research Training Grant ERB-FMBI-CT98-3438. 



Appendix 

We consider a group G with representation matrices t"^ = (t)j -^ . Then the quadratic Casimir 
C{R) of a representation R is defined by 

itH^)l = C{R)5l. (A.l) 

For SU{?>) triplets q and for SU{2)l doublets L we have 

Csui^{q) = \, Csu(2){L) = -^. (A.2) 

For U{1)y we have 

CU) = \y\f\ (A.3) 

where Y{f) is the hypercharge of the field /. The factor 3/5 is the grand unified normalisation. 
For the adjoint representation of the group of dimension d{G) we have 

C(G)5^^ = /^^^/^^^, (A.4) 
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where f^^'" are the structure constants. Specifically for the groups we investigate 

C{SU{3)c) = 3, C{SU{2)l) = 2, C(f/(l)y) = 0, (A.5) 

and C{SU{N)) = N. The Dynkin index is defined by 

Trfi(t^t^) = S{R)6^''. (A.6) 

For the respective fundamental representations / we obtain 

SU{3), SU{2) : S{f) = ^, (A.7) 

U{1)y : S{f) = h\f), (A.8) 

where we have inserted the GUT normalisation for f/(l)y. 

The coefficients in the two-loop running of the gauge couplings ( p.8|) are given by [^ 

B^:^ = (f,l,-3), (A.9) 
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